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We propose that in QCD with dynamical quarks, colour deconfinement occurs when an external field induced 
by the chiral condensate strongly aligns the Polyakov loop. This effect sets in at the chiral symmetry restoration 
temperature T x and thus makes deconfinement and chiral symmetry restoration coincide. The predicted singular 
behavior of Polyakov loop susceptibilities at T x is shown to be supported by finite temperature lattice calculations. 



1. INTRODUCTION 

Hot QCD matter undergoes two well-defined 
phase transitions at zero baryon chemical poten- 
tial depending on the bare quark mass parameter. 
For (m q = oo ) there is the deconfinement tran- 
sition associated with the spontaneous breaking 
of the global center symmetry Zn c 6 SU (N c ) of 
the pure SU(N C ) gauge theory. N c is the num- 
ber of colours of quarks. The order parameter 
describing this transition is the expectation value 
of the Polyakov loop L(T). There exists a critical 
temperature T^, with L(T) — V T < Td and 
L(T) > V T >T d . This transition is similar to 
the magnetic transition in Zn spin theories with 
the magnetization m(T) being analogous to L(T). 
The deconfinement transition and the transition 
in a spin system have same order when the num- 
ber of colours N c and the number of components 
of spin N are equal M . The nature of the decon- 
finement transition changes with N Cl for N c = 2 
the transition is second order and for N c > 3 it is 
first order. 

On the other hand for m q = 0, there is the 
chiral transition at the critical temperature T x 
associated with the restoration of chiral symme- 
try which is spontaneously broken at low tem- 
peratures. The order parameter describing this 
transition is the chiral condensate K = (tptp). 
The nature of the chiral transition depends on 
the number of quark flavours Nf. For Nf = 2 
the transition is second order and first order for 
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Nf > 3. In the case of a second order chiral tran- 
sition the critical behaviour of the chiral conden- 
sate are described the critical exponents defined 
in the following. 

K(T,m q = 0)~(T x ~Tf, T<T X (f) 

while at T — T x , it vanishes as 

K(T X — 0, m q ) ~ m q ^ s (2) 

when m q — > 0. The chiral quark mass suscepti- 
bility 

x£(T,m q ) = (^f) ~ ((#) a )-(V$ 2 

\ OTflq J rp 

~ |t-t x |-t (3) 

diverges at T = T x for m q — > 0. j3, S, 7, ... 
are the critical exponents. It is believed that for 
Nf — 2 the chiral transition is in the universality 
class of 0(4) according to symmetry arguments 
H , however it is yet to be confirmed by numerical 
experiments 

Both the symmetries of the Lagrangian, the 
global Zjv c symmetry and the chiral symmetry, 
are explicitly broken for finite non-zero m q . The 
role of m q for chiral transition is similar to that of 
an external magnetic field in spin system, so that 
for large m q , the temperature variation of the chi- 
ral condensate becomes completely smooth, like 
the magnetization m(T) in spin systems. 

It is not clear what role m q plays for the de- 
confinement transition. One might consider the 
inverse bare quark mass to play such a symmetry 
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breaking role for deconfinement, which would im- 
ply that for m q — > 0, the temperature variation of 
the Polyakov loop expectation value would then 
become smooth In lattice studies it is found, 
however, that this is not the case ||: even for 
m q — > 0, the Polyakov loop varies sharply with 
temperature and the corresponding susceptibility 
((L 2 ) — (I/) 2 ) peaks sharply at the chiral restora- 
tion temperature T x , which is considerably lower 
than Td- In somewhat loose terminology, one 
describes this situation by saying that in QCD 
for m q — > 0, deconfinement and chiral symmetry 
restoration coincide. 

To understand the underlying reason for this 
coincidence of the critical temperature of chiral 
and deconfinement transitions we propose Qthat 
the explicit breaking of the Zn symmetry is due 
to the effective external field which is inverse of 
the constituent quark mass ttiq |7| . It is believed 
that the constituent quark mass depends on the 
chiral condensate K, however, the exact depen- 
dence is not yet known. If the non-Goldstonc 
hadron masses are related to the chiral conden- 
sate, e.g., if the nucleon mass is given by ||] 

M n ~ 3 toq ~ K 1 ' 3 , (4) 

then chiral symmetry restoration with K — » 
gives rise to a sudden large increase of the exter- 
nal field H ~ 1/rriQ, forcing the Polyakov loop 
to take a large non-zero value. We want to argue 
that it is this effect which causes a sharp varia- 
tion of the Polyakov loop L(T) at T x < T d , with 
the temperature dependence of L(T) governed by 
that of the chiral condensate K(T), which is in 
general different from that obtained in pure gauge 
theory at Td- 

Although we believe that this interplay of chiral 
symmetry will in general do not depend on the 
value of Nf and N c , for simplicity and to compare 
with lattice results, we consider Nf — 2 and N c = 
3 ||. As we have already mentioned in this case 
the chiral transition is second order in the chiral 
limit, and deconfinement is first order in the limit 
of infinite heavy quarks. 

As indicated above, we now take the Zjq struc- 
ture of QCD to be that of Zjy spin theory in an 
external field H, determined by the constituent 
quark mass. From Eq. (|j), H and hence also 



the Polyakov loop thus become functions of K: 
the quark mass dependence of L enters through 
K(T,m q ), with L(T,K). In spin systems with 
external field, the magnetization m(T,H) is for 
H an analytic function of T and H . We there- 
fore assume that L(T, K) is for oo > K > an an- 
alytic function of T and K, with the quark mass 
dependence of L entering through K(T, m q ). 

Once the effective external field for deconfine- 
ment depends on the chiral condensate K, the 
critical behavior of the chiral condensate will be 
reflected in the behavior of the Polyakov loop at 
T = T x . For a first order chiral transition with 
a jump in the chiral condensate across the chiral 
transition, there will be a discontinuity in the ef- 
fective external field giving rise to a discontinuity 
in the Polyakov loop across the chiral transition 
point To go further one can write down 

for the variation of Polyakov loop L, 

WM dK. (5) 



<IL = \^-) dT 
\9Tj K 



dK 



From this we can calculate the variation of the 
Polyakov loop with respect to m q and tempera- 
ture T as follows, 

dL \ fdL\ f dK 
dK J rj, \ dm„ 



and 



dm 

dL_ 
df 



(G) 



di_ 

df 



K 



8L \ 
d~K J T 



dK 
~df 



■(7) 



Both these response functions of the Polyakov 
loop diverge at T = T X in the chiral limit m q = 0, 
since the chiral susceptibilities x£ = (dK /dm q )T 
and x£ = (dK /dT) mq diverge in this limit. 

In ||, the chiral susceptibilities and x^ 
were studied on an 8 3 x 4 lattice for quark masses 
m q a = 0.075, 0.0375 and 0.02. Results for x£ 
and x£ are shown in Figs. 1 and 2, respectively, 
as functions of the effective temperature variable 
T = k = 6/g 2 , where g denotes the coupling in 
the QCD Lagrangian. The peak heights of both 
these response functions increase with decreasing 
quark mass, which indicates the divergence in the 
chiral limit m q — ► 0. In Figs. 3 and 4, we show 
the corresponding results for the Polyakov loop 
susceptibilities. They also peak sharply, the 
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Figure 1. The chiral temperature susceptibility 
x K as function of the temperature variable n — 
6/5 2 - 



Figure 3. The Polyakov loop temperature suscep- 
tibility function of the temperature variable 
k = 6/g 2 . 
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Figure 2. The chiral quark mass susceptibility x£ 
as function of the temperature variable k = 6/g 2 - 



Figure 4. The Polyakov loop quark mass suscep- 
tibility x^ as function of the temperature variable 
k = 6/ 5 2 . 
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peak positions coincide with those for the chiral 
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susceptibilities, and here as well the peak height 



increases with decreasing quark mass. The ob- 
served behavior therefore provides clear support 
for the divergence of the temperature and quark 
mass derivatives of the Polyakov loop at T x in the 
chiral limit. 

A direct comparison of the functional behav- 
ior, ideally of the relevant critical exponents, for 
the two cases becomes difficult for two reasons. 
The chiral susceptibilities do not at present lead 
to the predicted 0(4) exponents; this may in- 
dicate that the quark masses are still too large 
for pure critical behavior. The Polyakov loop 
susceptibilities contain in addition unknown non- 
singular factors (dL/dT)K and (dL/dK)T, which 
will modify the non-singular m 9 -dependence rel- 
ative to that of the chiral susceptibilities. Lattice 
studies for smaller quark masses would therefore 
be very helpful. 

In summary, we have shown that chiral symme- 
try restoration, with the resulting sudden change 
in an effective constituent quark mass, leads to a 
suddenly increasing external field which aligns the 
Polyakov loops at T = T x < Td and thus produces 
a strong explicit breaking of the Zn symmetry of 
the gauge field part of the Lagrangian. This effect 
makes chiral symmetry restoration 'coincide' with 
deconfinement. The resulting predictions for di- 
verging Polyakov loop susceptibilities are found 
to be well supported by finite temperature lat- 
tice calculations for full QCD with two flavours 
of light quarks. 
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